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I.  INTRODUCTION 


We  consider  in  this  paper  the  derivation  of  governing  equations 
for  fluid  flows  through  ducts.  Such  flows  are  important  elements  in 
many  mechanical  systems.  Most  fluid  mechanics  textbooks  present, 
therefore,  a simple  derivation  of  the  governing  equations,  which 
reduces  the  general  three-dimensional  equations  to  a set  of  equations 
for  one-dimensional  flow.  Experience  has  shown  that  these  equations 
are  adequate  for  many  applications.  Probably  because  of  this  success 
researchers  sometimes  tend  to  disregard  the  limits  of  applicability 
of  the  one-dimensional  flow  equations.  In  order  to  derive  governing 
equations  for  more  complicated  flows  they  duplicate  the  steps  used  for 
simple  duct  flows.  The  resulting  equations  are  not  always  adequate, 
e.g.,  in  case  of  certain  non-steady  flows.  Some  textbooks  discuss 
limitations  of  the  usual  tube  flow  equations.  Often,  however,  the 
discussion  is  rather  general,  or  limited  to  examples  and  exercise 
problems,  and  easily  overlooked  by  casual  readers.  In  this  paper  we 
will  concentrate  on  the  limitations.  We  will  keep  the  discussions 
simple  by  considering  in  detail  only  a one-phase  flow  in  a straight 
duct  with  a constant  cross-section.  The  discussion  of  the  example  will 
provide  a methodical  approach  to  the  derivation  of  flow  equations  for 

^ more  general  cases. 

I 

The  starting  point  of  our  discussion  is  the  set  of  general  three- 
dimensional  flow  equations.  In  order  to  make  this  paper  self-contained, 
we  list  the  equations  in  Section  2.  In  Section  3 we  specialize  the 
equations  for  the  case  of  a duct  flow  using  a standard  procedure,  which 
is  found  in  textbooks.  In  order  to  establish  limits  for  the  validity 
of  the  specialized  equations,  we  carry  out  in  Section  4 a more  careful 
derivation  of  the  duct  flow  equations.  This  derivation  provides  quan- 
titative information  about  the  errors  which  are  introduced  by  the  special- 
ization of  the  equations.  A comparison  of  the  derivations  and  results 
of  Sections  3 and  4 reveals  that  in  standard  derivations  of  the  equations 
some  non-zero  terms  are  neglected.  In  Section  5 two  examples  are  pre- 
sented: a steady  flow  and  an  approximation  to  an  interior  ballistics 
flow.  Quantitative  estimates  are  given  for  some  usually  neglected  terms 
in  the  governing  equations.  Section  6 contains  some  conclusions  which 
can  be  dravm  from  the  discussions  o.*;  the  equiitions. 


2.  BASIC  GOVERNING  EQUATIONS 

We  consider  flows  which  satisfy  conservation  laws  for  mass,  momen- 
tum and  energy.  Governing  equations  for  such  flows  are  derived  and 


discussed,  e.g.,  by  Tsien  in  Reference  1 and  Batchelor  in  Reference  2. 

In  this  section  we  summarize  the  equations  in  order  to  ms’.o  this  paper 
self-contained.  We  use,  in  general,  the  same  notation  as  Tsien,  includ- 
ing the  convention  about  the  summation  over  equal  indexes. 

First  we  will  consider  the  equations  in  integral  form.  In  these 
equations  the  volume  integrals  are  for  an  arbitrary  control  volume  V, 
which  need  not  be  simply  connected.  We  assume,  however,  for  simplicity 
that  its  surface  S has  everywhere  an  outward  pointing  normal  nj . The 
conservation  of  mass  can  then  be  expressed  by  the  equation 


cl_ 

3t 


/ 


P 


dV  + 


)U.n .dS 
J 3 


0 


(2.1) 


The  momentum  equations  are 


1^  / pu^dV  + ^pu^u.n.dS  + / 1^  dV  = / Fj^dV.  (2.2) 

The  specific  kinetic  energy  of  the  fluid  is 


k = 


1 

2 


ujuy 


(2.3) 


Combining  eqs.  (2.1)  and  (2.2)  we  obtain  for  the  kinetic  energy  the 
equation 


3t 


J’ 0 k dV  i-^p  k u^n^dS-f^ 


i£_ 

|3x. 


dV 


= f U.F.; 

J 3 3 


(2.4) 


^H.S.  Tsien J "The  Equations  of  Gas  Dynamics y"  in  Fundamentals  of  Gas 
DynamioSy  edited  by  H.W.  Emmons,  Princeton  University  Press,  1958. 

^G.K.  Batchelor,  An  Introduction  to  Fluid  Dynamics,  Cambridge  University 
Press,  1967. 
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The  first  law  of  thermodynamics  is 


3t 


y*  p e dV+^( 


e 


(2.5) 


By  adding  eqs.  (2.4)  and  (2.5)  we  obtain  an  equation  for  the  specific 
total  internal  energy  e + k: 


3t 


J Pie  + k)dV+^ip 


(e  +.k)u.n.dS  + 
J .1 


u^n^dS  =J^-|^^  dV +y((|)+  UjFj)dV. 


(2.6) 


The  last  integral  in  eq.  (2.6)  is  the  contribution  of  external  and 
viscous  forces  to  the  changes  of  the  total  internal  energy.  Its  first 
part,  /<|)dV,  is  the  contribution  of  viscous  forces  to  the  internal  energy 
e.  The  integrand  $ , i.e.,  the  heat  dissipation  function,  can  be  ex- 
pressed in  terms  of  the  viscous  stress  tensor  Tjic: 


$ = 


(2.7) 


We  assume  that  the 
mula  to  the  strain  rate 
page  132) 


viscous  stress  tensor  is 
tensor  ej^j  (Reference  1, 


related  by  Stokes  for- 
page  13,  Reference  3, 


kj 


= 2u  G 


(f'  - 


2 


6 , .e . . 
kj  11 


(2.8) 


This  definition  is  not  restricted  to  constant  viscosities  p and  p',  i.e., 
to  homogenous  fluids.  However,  it  restricts  the  considerations  to  iso- 
tropic fluids.  The  viscosities  p and  p*  must  be  positive  or  zero. 


Hamel,  Meohanik  dev  Kontinm,  R,G.  Teubner,  Stuttgart,  1956. 
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The  strain  rate  tensor  is  defined  by 


(Reference  2,  page  80) 


^kj 


2 ^3x. 
3 


(2.9) 


Substituting  (2.8)  and  (2.9)  into  eq.  (2.7)  we  obtain  the  following 
expressions  for 


"kj'kj  * 3"^  \j  = ii=kj 


= 2\i  {e^} 


trace 


(u* 


2 1 


{c}2 


trace 


3u,  3u.  2 

1 K 


(2.10) 


Eq.  (2.10)  shows  that  the  heat  dissipation  function  ^is  always  positive 
for  a stress  tensor  of  the  form  (2.8). 

The  term  /ujF-dV  in  eq.  (2.6)  is  the  contribution  of  viscous  and 
body  forces  to  tne"^ changes  of  the  kinetic  energy  k.  The  force  (per  unit 
volume)  Fj  is  a sum  of  body  forces  pXj  and  viscous  forces  Tj . The  latter 
can  be  expressed  in  terms  of  the  viscous  stress  tensor  We  thus 

have  the  equation  ^ 


F. 

1 


= pX.  + T. 
3 3 


Combining  eqs.  (2.7)  and  (2.11)  we  obtain 


$ + u.F.  = 6 + u.T.  + pu.X.  = 
^ 3 3 ^ 3 3 3 3 


(2.11) 


(2.12) 
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In  this  form  we  have  subdivided  the  contributions  of  forces  to  the  ' 
changes  of  the  total  internal  energy  into  contributions  by  viscous  and 
by  body  forces.  The  corresponding  volume  integral  in  eq.  (2.6)  is 


u.FJdV  = f($+  u.T.  + pu.X.)  dV  = 
1 1 J ^ 3 3 3 3 

= (TT-.u.n,  dS  + Ipu.X.dV 
J 3k  3 k ^ J 3 3 


(2.13) 


We  note  that  the  surface  integral  in  eq.  (2.13)  contributes  to  the  inter- 
nal as  well  as  to  the  kinetic  energy  of  the  flow.  It  represents  the 
viscous  forces  acting  on  the  surface  of  the  control  volume.  The  volume 
integral  over  the  body  forces  contributes  to  the  kinetic  energy  only. 


The  governing  equations  (2.1),  (2.2)  and  (2.4)  through  (2.6)  can 
also  be  expressed  in  differential  form  as  follows 


3 


3x. 

3 


k ’ 


Ij-  (p  k)  » ^ Cp  k u.)*  u.|L=  u.F.  . 


3 


3u. 


3t 


(P  0)  < (P  e u ) . p Q - 

3 3 3 . 

a a 

r?(e+k)i  + ^ [p(e+k)Uj]  + (P  Uj)  = 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 


Eq.  (2.16)  is  a consequence  of  eqs.  (2.14)  and  (2.15),  because  k is 
defined  by  eq.  (2.3).  Also,  eq.  (2.18)  is  the  sum  of  eqs.  (2.16)  and 
(2.17).  We  have,  therefore  only  five  independent  differential  equa- 
tions for  the  six  quantities  p,  uj^,  p,  and  e.  To  complete  the  system 
of  equations  we  need  another  equation,  which  is  provided  by  the  equa- 
tions of  state  for  the  fluid  under  consideration.  We  assume  that  such 
an  equation  is  available,  e.g.,  in  the  form 


0(e,p,p)  = 0 , 


(2.19) 
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and  that  eq.  (2.19)  can  be  solved  explicitly  for  either  of  the  argu 
ments.  For  example,  in  case  of  an  ideal  gas  with  constant  specific 
heats  eq.  (2.18)  is 


JL-  (Y-l)e  = 0 . (2.20) 

P 


For  the  discussions  in  the  rest  of  this  paper  we  will  not  make  use  of 
eq.  (2.19)  or  (2.20).  The  assumption  of  the  existence  of  such  an  equa- 
tion is  made  here  only  to  close  the  set  of  governing  equations. 


3.  APPROXIMATE  GOVERNING  EQUATIONS  FOR  DUCT  FLOWS 

In  this  section  we  derive  approximate  governing  equations  for  duct 
flows.  The  dominant  component  of  such  flows  is  usually  in  the  axial 
dire^'tion.  Also,  in  many  cases  only  the  dependence  of  flow  properties 
on  the  axial  coordinate  is  of  practical  interest.  Duct  flows  are  there- 
fore usually  treated  by  one-dimensional  equations  which  are  derived 
from  the  general  flow  equations  of  Section  2. 

A standard  procedure  for  the  iei'ivation  of  these  equations  is  to 
consider  a control  volume  which  consists  of  a length  Az  of  the  duct. 

The  integral  forms  of  the  governing  equations  are  applied  to  this 
control  volume  and  corresponding  differential  equations  obtained  by 
letting  Az  approach  zero.  This  method  is  used,  e.g.,  in  References  4 
and  S,  and  we  will  follow  these  references  closely. 

Another  possible  approach  is  to  start  with  governing  equations  for 
one-dimensional  flow,  i.e.,  a flow  which  depends  on  only  one  coordinate 
and  which  has  a velocity  component  in  the  direction  of  that  coordinate 
only.  Three  dimensional  effects,  e.g.,  from  the  wall  friction,  are  then 
added  to  the  equations  by  ad  hoc  procedures.  We  will  not  pursue  this 
approach  here  because  the  former  approach  can  be  generalized  more 
easily. 


^A.H.  Shapiro,  The  Dynamiae  and  ThermodynamioB  of  CompreBsible  Fluid 
Flow,  Vol.  I and  II,  Roland  Press  Company,  New  York,  1954. 

^L.  Crooao,  "One-Dimensional  Treatment  of  Steady  Gas  Dynamos"  in 
Fundamentals  of  Gas  Dunamias  II.  edited  by  H.N.  Emmons,  Princeton 
University  Press,  1958. 
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Let  z be  the  axial  coordinate  and  let  for  simplicity  the  cross- 
sectional  area  A of  the  duct  be  constant.  The  continuity  equation  (2.1) 
is  then  for  the  control  volume 


The  bars  on  P and  u in  eq.  (3.1)  indicate  that  we  are  dealing  with 
average  density  and  velocity,  respectively.  We  apply  now  the  mean  value 
theorem  to  the  first  term  in  eq.  (3.1)  and  use  a Taylor  series  expansion 
for  the  second  term.  The  result  is 


{p(z)  A}  • Az  + Az  • A = 0(Az2)  , (3.2) 


where  z < z < z + Az.  Letting  Az  in  eq.  (3.2)  approach  zero  we  obtain  the 
continuity  equation 


3t  3z  " ^ • 


(3.3) 


The  momentum  balance  equation  is  considered  in  the  z-direction  only. 
First  we  obtain  as  above  from  eq.  (2.2) 


_ _2 

{p(z)  u(z)  A}‘Az  + Az*A  = F*A‘Az  + O(Az^) 


(3.4) 


The  momentum  equation  for  the  average  flow  properties  is  obtained  from 
eq.  (3.4)  by  letting  Az  approach  zero.  The  result  is 


^ (p  D)  + ^ (p  u^) 


(3.5) 
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The  force  per  unit  volume,  F,  can  be  expressed  as  a sum  of  two  compo- 
nents in  analogy  to  eq.  (2.11).  The  momentum  equation  is  then 


3 3 _2, 

^ (P  a)  + ^ (P  u ) 


p X + T . 


(3.6) 


The  quantity  T in  eq.  (5.6)  is  obtained  from  the  resultant  of  the  viscous 
boundary  forces  on  the  surface  of  the  control  volume.  For  simple  tubes 
T can  be  expressed  in  terms  of  the  pipe  friction  coefficient  or,  by  exper- 
imental correlations^  iji  terms  of  the  surface  roughness  and  perimeter  of 
the  tube.  The  term  p X usually  represents  the  gravity  force  component  in 
the  axial  direction  of  the  tube. 

A combination  of  eqs.  (3.3)  and  (3.6)  yields 


^ (j  p u^)  + --  (j  p u^)  + u II  = u F = u p X + u T . (3.7) 


Eq.  (3.7)  can  be  considered  as  an  equation  for  the  kinetic  energy,  if 
the  latter  is  approximated  by 


k = 


(3.8) 


However,  eq.  (3.7)  is  a mathematical  consequence  of  the  continuity  and 
momentum  equations,  i.e.,  eqs.  (3.3)  and  (3.6),  and  is  independent  of 
any  assumptions  about  the  kinetic  energy. 

Next  we  consider  the  energy  balance.  Following  general  practice 
(see,  e.g..  Reference  4)  we  start  with  the  eq.  (2.6)  instead  of  using 
the  first  law  of  thermodynamics,  i.e.,  eq.  (2.5).  For  the  control  vol- 
ume we  obtain  first 

{p*(e+lc)*A)*Az  + AZ'A  -r-  {p*u*(e+T()  + u«p}  = 

ot  0 Z 

= (Q  - |^)*A‘Az  + ($+  u*F)*A*Az  + O(Az^)  . (3.9) 
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At  the  limit  Az  -+■  0 eq.  (3.9)  yields  the  energy  equation 

{p  (e+ic)}  + {p  u (e+ic)  + u p}  = 

= + $ + (3.10) 

oZ 

i'he  equation  of  state,  such  as  eq.  (2.19),  cont£ins_usually  the 
internal  energy  e and  not  the  total  ijiternal  energy  e + k.  Therefore, 
eq.  (3.10)  is  modified  to  eliminate  k.  To  this  end  it  is  assumed  that 
the  approximation  (3.8)  holds,  and  eq.  (3.7)  is  subtracted  from  eq. 

(3.10) . The  result  is 


(P  e)  + (P  u e)  + p = Q - II  ^ . (3.11) 


Eq.  (3.11)  is,  of  course,  the  first  law  of  thermodynamics  and  could  have 
been  obtained  directly  from  eq.  (2.5)  without  any  assumptions  about  the 
kinetic  energy. 

In  order  to  use  eqs.  (3.3),  (3.7),  and  (3.n^  for  imputations  we 
need  among  other  data^estimates  for  the  forces  p X and  T and  for  th£  heat 
dissipation  function  The  latter  is  often  expressed  in  terms  of  T by 
the  following  arguments.  (See,  e.g..  Reference  4,  page  39  ff.  and 
972  ff.) 

The  last  integral  on  the  right  hand  side  of  the  energy  equation 
(2.6)  is  according  to  eq.  (2.13) 


W = / (*+  u.F.)dV  = ^T..u.n.dS  + / p u.X.dV  . (3.12) 


The  surface  integral  in  eq.  (3.12)  represents  the  work  rate  of  viscous 
forces  acting  on  the  surface  S of  the  control  volume.  We  subdivide  this 
surface  into  material  boundaries  (e.g.,  duct  walls)  and  flow-through 
surfaces  Sq.  The  work  done  on  material  boundaries  is  called  shaft  work. 
The  work  by  viscous  forces  on  the  flow-through  boundaries  is  called 
shear  work.  Let  the  corresponding  work  rates  be  Wg^aft  Wq,  respec- 
tively. In  these  terms  eq.  (3.12)  is 


W = / ((|)  + = W 


shaft 


"o  * / P 


U.X.dV  . 


(3.13) 
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The  integral  Wq  over  the  open  boundaries  is  usually  assumed  to  be  negli- 
gible. For  example,  in  case  of  a tube  flow  it  is  argued  that  integration 
over  the  core  flow  region  contributes  very  little  to  the  integral  because 
is  small  in  that  region.  Integration  over  the  boundary  layer  region 
also  contributes  little  because  the  velocity  Uj  is  small  in  the  boundary 
layer.  Hence  Wq  must  be  small. 

If  we  carry  out  the  derivation  of  eq.  (3.10)  using  the  relation 
(3.13)  we  obtain 

l^tP  C5*k)!  • |j(5  D (rt)  . u p)  = 5 - 11 . . 5 a )!  . (3.14) 


Combining  eqs.  (3.14),  (3.7),  and  (3.8)  we  obtain  as  the  first  law  of 
thermodynamics  instead  of  eq.  (3.11)  the  equation 


Cp  5)  * C5  u 5)  . T II  = Q - II  . - u f , (3.15) 

In  this  equation  the  heat  dissipation  function  $ is  approximated  by 


= + W -tiT. 

shaft  0 


(3.16) 


At  the  material  boundaries  the  velocity  of  the  fluid  is  equal  to  the 
velocity  of  the  boundary.  Therefore,  Wshaft  is  non-zero  only  if  the 
boundaries  are  moving.  If  the  tube  does  not  contain  moving  boundaries 
and  ¥q  is  neglected,  then  eq.  (3.16)  becomes 


$ = -n  T, 


A / J 


which  is  the  usual  approximation  of  $ for  tube  flows  (Reference  4, 
page  972  ff.) 


_ In  case  of  two-phase  flows,  e.g.  particles  submerged  into  the  fluid, 
Wshaft  i^  assumed  M be  the  work  of  drag  forces.  I^et  the  average  par- 
ticle velocity  be  Uparticle  ^he  drag  force  be  Then 


shaft 


TI  T 

particle  drag  ’ 


(3.18) 
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The  resultant  T of  the  viscous  forces  is  in  this  case  the  sum  of  par- 
ticle drag  and  wall  friction  forces 


T = 


wall 


(3.19) 


The  equation  for  becomes  then 


$ ” ^^particle  ^drag  ” ^ ^wall  ' 


This  equation  is  sometimes  modified  by  an  ad  hoc  factor,  see  Reference  6, 
page  81. 

In  summary,  either  eq.  (3.17),  or  eq.  (3.20)  provides  a convenient 
estimate  for  $ i.e.,  for  the  rigfit-hand  side  of  the  energy  equation 
(3.15).  Estimates  of  and  T^j  are  also  needed  to  express  the 

terms  on  the  right-hand  sides  of  tne^momentum  equation  (3.6).  It  appears 
from  the  derivation  that  no  further  estimates  of  flow  properties  are 
needed  under  quite  general  conditions. 

Some  limitations  of  the  approximation  (3.17)  become  obvious  if  we 
consider  non-steady  fluctuating  flows.  In  such  flows  it  is  possible 
that  the  signs  of  T and  u are  temporarily  equal.  In  these  cases  Wq 
cannot  be  neglected,  because  otherwise  we  would  have  a negative  heat 
dissipation  function.  Thus  it  seems  appropriate  to  ask  how  accurate  is 
the  energy  equation  (3.15).  Our  derivation  does  not  provide  any  clues 
to  an  answer  to  this  question.  We  will  therefore  rederive  the  duct  flow 
equations  more  carefully  in  the  next  section,  keeping  track  of  all  approx- 
imations involved. 


4.  PRECISE  GOVERNING  EQUATIONS  FOR  DUCT  FLOWS 

In  this  section  we  will  derive  complete  one-dimensional  governing 
equations  for  flows  through  constant  area  ducts,  including  formulas  for 
quantities  which  were  neglected  in  Section  3.  We  will  then  discuss  the 
differences  between  the  complete  equations  and  those  of  the  previous 


^G.B.  WalliBj  One-Dimenaioml  Two-Phase  Flow,  MoGvaw-Hill  Co,  New  York, 
1969. 
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section,  indicating  where  the  previous  derivation  of  the  equations  is 
insufficient. 


One-dimensional  duct  flow  equations  are  relations  between  average 
flow  properties.  The  equations  depend  therefore,  among  other  things,  on 
the  definitions  of  the  averages.  For  steady  duct  flows  certain  averages 
and  corresponding  governing  equations  have  been  discussed  by  Crocco  in 
Reference  5.  Because  the  averages  defined  by  Crocco  cannot  be  used  for 
non-steady  flows,  our  analysis  will  be  different.  The  results  of  our 
analysis  can  be  applied  to  steady  as  well  as  non-steady  flows. 


First  we  consider  the  continuity  equation  (2.1). 
volume  which  consists  of  a length  of  of  the  duct. 


For  a control 
eq.  (2.1)  is 


V ^ . 

" * 

[ jo  * 

jp  Ujds 

X A 

Am 

. A 

X3.AX3 


= 0 


(4.1) 


The  integrals  /pds  and  Jpu-^ds  are  functions  of  X3.  We  expand  these 
functions  in  Taylor  series,  intechange  the  order  of  integration  over  X3 
and  differentiation  with  respect  to  t in  the  first  term  of  eq.  (4.1), 
and  apply  the  mean  value  theorem  to  that  term.  The  result  is 


‘>^3  [If  i " 

dsl 

J x,+  6Ax- 

* ‘’'3  [ir  i ^ “3“’=] 

(4.2) 

A 

with  0 < 0 < 1. 

3 3 

A X3 

At  the  limit 

Axj  -»■  0 eq. 

(4.2)  yields 

3 f 

“3“= ' “ • 

(4.3) 

A ^3  A 


We  now  define  for  each  cross-section  X3  = const,  an  average  fluid 
density  p by 


(4.4) 


and  an  average  fluid  velocity  u by 


(4.5) 


The  continuity  eq.  (4.3)  can  then  be  expressed  in  terms  of  the  average 
density  and  velocity  as 


3p 

3t 


+ 


0 


(4.6) 


Eq.  (4.6)  is  identical  to  the  continuity  equation  (3.3).  However, 
we  have  now  established  that  th£  continuity  equation  is  of  this  form 
only  if  the  average  quantities  p and  u are  defined  by  eqs.  (4.4)  and 
(4.5),  respectively.  Thus,  if  we  chose  an  alternate  definition  of  the 
average  velocity,  e.g.,  the  simple  spatial  average 

'Jf=^/u.ds,  (4.7) 

A ^ J 


then  the  corresponding  continuity  equation  would  be 


3p  , 9 

3t  . 8X3 


(H  - u)]  . 


(4.8) 


The  right-hand  side  of  eq.  (4.8)  is  non-zero  in  general. 

The  momentum  equation  (2.2)  yields  for  a duct  flow  in  analogy  to 
eq.  (4.3) 


9 r 

A 


+ aT-  / P Uj^U3ds  + / ds  = / F^ds 


3x, 


(4.9) 


If  the  flow  is  axially  symmetric,  then  for  k = 1 and  k = 2 all  terms 
in  eq.  (4.9)  are  identically  zero.  In  cases  of  non-axisymmetric  flows 
all  three  momentum  equations  are  needed  to  describe  the  flow,  e.g.,  in 
the  case  of  a non-vertical  tube  in  a gravity  field.  We  will  consider 


for  simplicity  only  the  momentum  equation  in  the  axial  x^-direction, 
thus  restricting  the  analysis  to  axisymmetric  flows.  Eq,  (4.9)  thus 
becomes 


Ujds  . 


_3 

3x 


/ P 
3 A 


u,ds 

O 


I = / 

A ^3  A 


F3ds 


(4.10) 


Eq.  (4.10)  contains  two  new  flow  variables  for  which  averages  have  to 
be  defined.  We  chose  the  following  definitions: 

p = ^ /p  ds  (4.11) 

^ A 

and 


Expressing  the  momentum  equation  (4.10)  in  terms  of  average  quantities 
we  then  obtain 


with 


I^Cp  -u)  . ^(-p  -u^)  . f- . (4.13) 

= il:  ' 

A 

(4.14) 

“ ^3  { “ “HI 


Comparing  the  momentum  equation  (4.13)  with  the  corresponding 
eq.  (3.5),  we  see  that  the  latter  equation  is  in  error.  The  reason  for 
this  error  is  that  eq.  (3.4)  should  have  contained  the  term  Az.Cm'A.  Eq. 
(4.14)  shows  that  this  term  is  non-zero  in  general.  The  expression  in 
square  brackets  in  eq.  (4.14)  is  negative  or  zero  according  to  Schwarz's 
inequality.  It  is  zero  if  and  only  if  U3  = const,  across  the  duct. 

Hence  the  correction  term  Cj^  is  zero  only  in  case  of  a slug  flow  or  if 
the  term  is  independent  of  X3.  The  latter  is  the  case  for  steady  incom- 
pressible flows  through  constant  area  ducts.  In  all  more  interesting 
cases  Cm  is  non-zero  and  its  magnitude  should  be  estimated  to  justify 
the  neglect  of  C^,  or  C^  should  be  included  in  the  momentum  equation. 
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The  average  force  per  unit  volume,  F3,  is  defined  by  eq.  (4.12). 
For  later  reference  we  note  that  according  to  equation  (2.11)  F3  is  a 
sum  of  body  forces  and  viscous  forces.  We  define  the  corresponding 
averages  by 


and 


(4.15) 


With  these  definitions  eq.  (4.18)  becomes 


(4.22) 


where 


r - ri 

el  P 3x, 


■M 


(4.23) 


u)  - \ 


e u^ds 


(4.24) 


The  nature  of  the  correction  terms  Cg^  and  Ce2  is  similar  to  that  of  the 
correction  term  in  the  momentum  equation.  They  are  zero  for  slug  flow 
and  should  be  estimated  in  other  cases.  If  we  compare  eq.  (4.22)  with 
the  corresponding  eq.  (3.11),  we  see  that  the  latter  is  in  error.  The 
reason  for  the  error  is  an  oversight  of  a term  4z(Cg2  + Cg2)*  A which 
should  have  been  introduced  in  eq.  (3.9).  The  correction  terms  enter 
the  equations  because  a product  of  function  averages  is  in  general  not 
equal  to  the  average  of  the  product  of  the  functions.  Or,  differently 
expressed,  multiplications  of  functions  and  averaging  of  functions  are 
not  commutative  operations. 

We  mentioned  in  Section  3 that  (|>  is  usually  approximated  by  -u  T3. 

It  was  also  shown  that  such  an  approximation  is  based  on  the  assumption 
that  a term  Wq  can  be  neglected.  We  will  now  investigate  the  approxima- 
tion more  carefully.  By  the  definition  (4.20)  we  have 


■kj  3X, 


A \ 3x.  ■ A j ^ 3x^  “ 

A A 


(4.25) 


The  first  term  on  the  right-hand  side  of  eq.  (4.?.3)  we  recognize  as 
Wq,  i.e.,  the  average  of  the  gradient  of  the  work  rate  of  viscous 
forces_on  the  cross  section  A.  The  second  term  may  be  approximated 
by  -u  T3.  The  final  formula  for  ^ is  then 


where 


$ = - u Tj  ♦ Wq  . , 

3 A 


and 


C 


(4.26) 


(4.27) 


(4.28) 


Combining  eqs.  (4.26)  and  (4.22),  we  obtain  for  the  energy  equation 
(first  law  of  thermodynamics)  the  expression 


3 

at 


(P  e) 


a 

9x, 


_ an 

(P  e u)  * p ^ 


H - 


u T + C , 
3 el 


+ C ^ + W + C. 
e2  0 4 


(4.29) 


The  first  two  correction  terms,  Cg^  and  Co*  appear  in  eq.  (4.29) 
because  of  the  averaging  of  some  terms.  Tn£  las'^^o  correction  terms. 

Wo  and  0$,  are  due  to  the  approximation  of  $ by -u  T3. 

The  equation  for  the  kinetic  energy  can  be  treated  formally  in 
the  same  manner  as  the  equation  for  the  internal  energy.  One  can  intro- 
duce error  terms,  corresponding  to  Cg^  and  Cg2>  either  in  the  kinetic 
energy  equation  or  in  the  equation  for  the  total  internal  energy. 

Since  typically  only  the  equation  for  internal  energy  is  needed  for 
computations,  the  other  equations  are  not  formally  derived. 

In  summary,  we  have  shown  that  the  one-dimensional  governing 
equations  for  average  flow  properties  in  duct  flows  are  not  the  same  as 
equations  for  locally  one-dimensional  flows.  If  the  medium  is  compressible 
then  the  additional  terms  in  the  governing  equations  vanish  only  for  slug 
flow.  For  other  flows  the  magnitudes  of  the  terms  should  be  estimated  for 
each  case  to  check  their  significance.  Formulas  given  in  this  section  may 
be  used  for  that  purpose. 
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If  the  duct  is  axially  symmetric,  it  is  more  convenient  to  use 
cylindrical  coordinates  than  the  cartesian  coordinates  of  this  section. 
We  give  therefore  in  Appendix  A all  pertinent  formulas  in  cylindrical 
coordinates . 


5.  EXAMPLES  OF  TUBE  FLOWS 

5.1  Incompressible  Steady  Flow  Through  Cylindrical  Tubes 

In  the  case  of  an  incompressible  steady  duct  flow  the  flow  velocity 
is  constant  along  the  duct  and  dependent  on  the  radial  coordinate  r only. 
Also,  only  the  axial  coordinate  u of  the  velocity  is  non- zero.  There- 
fore, of  all  the  correction  terms  given  in  Appendix  A,  only  C$  can  be 
non-zero  in  this  case.  It  is  given  by  eq.  (A. 35),  which  reduces  to 

R 

C.  = uT  - (rt  )dr.  (5.1.1) 

<!>  z J 3r  '•  rz"^ 

0 

The  shearing  stress  Tj.2(r)  is  in  the  present  case  a linear  function 
of  r.  This  is  a conseq^ience  of  the  second  momentum  equation  (A. 13) 
which  reduces  to 


The  left-hand  side  of  e.,. 
linear  in  r: 


iE. 11_ 

3z  ‘ r 3r 

(.5.1.2)  is 


(r  T ) . 
rz^ 


constant. 


(5.1.2) 

Therefore,  Tr2  must  be 


(5.1.3) 


Substituting  eq.  (5.1.3)  into  eq.  (5.1.1),  we  obtain 


C.  = u T - 
V z 


4t  (R) 
rz^ 


r dr  = 


" "z  - 2 " ^rz('^5  i • 


(5.1.4) 
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The  average  shear  stress  is,  according  to  eqs.  (A. 17)  and  (A. 21), 


T 

z 


(5.1.5) 


Substituting  eq.  (5.1.5)  into  eq.  (5.1.4)  we  see  that  the  correction 
term  €$  is  zero. 

Hence  the  average  flow  equations  are  exact  for  incompressible 
steady  flows  through  circular  tubes.  This  is  essentially  a consequence 
of  eq.  (5.1.3)  and  the  result  is  valid  for  either  turbulent  or  laminar 
flows.  Also,  we  have  not  made  use  of  Stokes  equations  for  the  stress 
tensor,  nor  made  any  assumptions  about  the  viscosity  of  the  fluid. 

5.2  Lagrange's  Interior  Ballistics  Flow 

As  an  example  for  non-steady  tube  flows  we  consider  Lagrange's 
approximation  to  interior  ballistics  flow  (Reference  7) . The  approxi- 
mation is  obtained  by  postulating  that  the  average  axial  velocity  u 
of  the  gas  in  a gun  tube  is  at  any  time  a linear  function  of  the  axial 
distance  z,  i.e.. 


where  Zp(t)  and  il  = dZp/dt  are  the  location  and  velocity  of  the  pro- 
jectile, respectively.  ^^’We  assume  that  the  local  velocity  can  have  axial 
as  well  as  radial  components  which  may  depend  on  z,  t,  and  on  the  radial 
coordinate  r. 

Some  consequences  of  the  assumption  (5.2.1)  are  discussed  in  the 
Appendix  B.  In  summary,  the  discussion  shows  that  this  assumption, 
complemented  with  a second  Lagrange's  assumption 


P 


P 


^t°) 

0 2p(t) 


(5.2.2) 


^J.  Corner^  Theory  of  the  Interior  BallietiQB  of  Guns,  John  Wiley  and 
SonSj  New  Yorkt  1950. 
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is  oonsisttiut  with  tiio  isverage  continuity  equation  (4.6)  for  the  flow. 
One  can  also  assume  tliat  tlio  local  velocity  vector  has  the  form 


u* 


u(z,t)‘f(r) 

v(z,t)‘h(r) 


(5.2.3) 


l-or  any  roasonal)lo  functions  u(z,t)  and  f(r)  one  can  determine  corres- 
ponding funtions  v(z,t)  and  h(r)  such  that  the  local  continuity  equa- 
tion is  satisfied.  (The  necessary  formulas  are  jjiven  in  Appendix  B.) 
However,  a flow  characterized  by  eqs.  (5.2.1)  through  (5.2.3)  in  general 
does  not  satisfy  the  local  momentum  equations  if  constant  viscosities 
are  assumed.  Hence  !,agrange's  approximation,  (5.2.1)  and  (5.2.2),  and 
a local  velocity  field  of  the  type  (5.2.3)  can  be  consistent  only  for 
inliomogeneous  media,  i.e.,  media  with  variable  viscosity. 

Because  an  exact  solution  of  the  viscous  tube  flow  equations  is 
not  available,  we  cannot  obtain  exact  values  for  all  correction  terms. 
However,  the  correction  term  in  the  momentum  equation  is  independent 
of  the  viscosities  and  can  be  computed  exactly  for  any  flow  profile. 

In  contrast,  the  correction  terms  in  the  energy  equation  can  be  computed 
only  if  additional  information  is  available  about  the  stress  tensor  and 
the  internal  energy  profile.  These  terms  we  will  estimate  by  computing 
tlu'ir  values  for  constant  viscosities  and  for  a number  of  "reasonable" 
flow  profiles.  We  expect  by  such  calculations  to  ol)tain  at  least  order- 
of-magni tilde  estimates  of  the  correction  terms. 

Particularly  wo  will  consider  flow  profiles  of  two  types,  first 
wo  will  assume  a flow  field  which  is  described  by 


and 


u(r,z,t) 


u (t)  - 

p^  ^ n-t-2 

z (t)  n 
P 


ujt) 


v(r,  = ,t)  = - 


/r\  " 

1 - 

(^)  J 

r , 

\ ^ 

1 1 

\r) 

L 

) (5.2.4) 


This  flow  field  has  a Hagen-Poiseuille  profile  for  n = 2.  For  larger 
values  of  n it  approximates  turbulent  flow  profiles  or  profiles  with 
thin  boundary  layers. 

As  a second  example  we  will  consider  a flow  profile  which  approx- 
imates the  universal  profile  for  steady  turbulen  tube  flow. 
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The  flow  field  defined  by  eq.  (S.2.4)  satisfies  the  local  continu- 
ity equation,  if  the  density  is  given  by  eq.  (5.2.2).  We  note  the  inter- 
esting fact  that  local  continuity  requires  the  radial  flew  component  to 
be  directed  toward  the  center  of  the  tube.  This  is  due  to  the  higher 
mass  flow  rate  at  the  center  and  due  to  the  assumed  increase  of  the 
average  axial  velocity  u(z,t)  with  z. 


The  correction  term  C^,  of  the  momentum  equation  is  given  for  our 
flow  by  eq.  (B.48) 


m 


Substituting 


(5.2.5) 


cf  ^ n+2 

f(r)  = -- 


* -'I 


(5,2.6) 


into  eq.  (5.2.5),  we  obtain 


C 

m 


1 ^ 
n+1  9z 


(p  u2). 


(5.2.7) 


The  momentum  equation  is 


therefore  in  terms  of  the  average  axial  velocity 


3(P  G)  , 
at 


i-  (p  02)  + _ i> 

az  az  z 


(5.2.8) 


Eq.  (5.2.8)  shows  that  in  the  case  of  a Hagen-Poiseuille  profile  the 
momentum  transport  term  in  the  momentum  equation  should  be  increased  by 
about  33%.  Even  for  a rather  flat  profile  with,  say,  n = 10  the  correc- 
tion term  is  9%  in  this  example. 

The  first  correction  term  Cgj  of  the  energy  equation  is  zero  in 
our  example  because  the  density  p is  independent  of  r and  z.  (Soe 
Appendix  B for  a discussion  of  this  term.) 
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The  second  correction  term  Ce2  the  energy  equation  is  (see  eq. 
(B.54)) 


c = — 
e2  3z 


(1-f)  e r dr 


p 2 3 


z 2 3z 
n R 


^ ^z  \ (1-f)  e r dr)  . 


C 0 


(5.2.9) 


In  order  to  compute  this  term  we  would  need  to  make  an  assumption  about 
the  internal  energy  function  e.  For  the  present  discussion  we  will  not 
make  any  assumptions  and  leave  eq.  (5.2.9)  unchanged. 

The  average  heat  dissipation  function  which  appears  on  the  right 
hand-side  of  the  energy  equation  can  be  computed  by  the  formulas  (B.56) 
and  (B.57).  The  result  of  the  computation  is 
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i-iil.  fir  (:  113.2 
1,2  n \i^j  t.*l  I J 


(5.2.10) 


The  liquation  for  the  average  internal  energy  (first  law  of  thermo- 
dynamics) is  in  our  case 


3t  3z 


(p  e u) 


+ P 


3u 

3z 


= H + $ + Cg2  • 


(5.2.11) 


Substituting  eqs.  (5.2.9)  and  (5.2.10)  into  eq.  (5.2.11),  we  obtain 

R 


9(P  g)  ^ 1_ 

at  32 


= H + 2ii  n2 


/“ 

p j (1-f)  e r d^ 


(n^2)^  + /i  nil  + 


_ 3u 
" P 37 


2n 


3 n+1  2mAz 


WJ 


(5.2.12) 


In  eq.  (5.2.12)  we  have  included  the  correction  term  Cg2  into  the 
energy  flux  terra  on  the  left  hand  side.  It  is  readily  apparent  from  the 
form  of  the  term  that  the  correction  is  zero,  if  the  specific  internal 
energy  e is  independent  of  r. 
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In  the  heat  dissipation  function  on  the  right-hand  side  of  eq. 

(5.2.12)  the  terra  with  the  factor  fR/z)^'  can  generally  be  neglected, 
because  (R/z)^  is  of  the  order  10-^.  (R/z  is  large  in  the  vicinity  of  the 
breech,  where  the  one-diraensional  approxiraation  should  not  be  used  any- 
wayO  The  other  term  in  the  square  brackets  is  usually  replaced  by 
-u  Tj.  If  this  is  done,  then  two  additional  correction  terms  should  be 
included  in  the  equation.  The  general  formulas  for  these  terms  are 
given  by  eqs.  (B.61)  and  (B.62).  They  are  in  our  case 


and 


-ftj  [?  i ■ • ©'  ■■  • 


r dr  + y'  - j P) 


. lo- 


n+1  ■ 3 


•] 


'u'\2 

Py 


C.  = ut^»(j£1  1. 


„2^,2  (RV  v2  2., 2 

R h'  + -r  h + z f' 


(5.2.13) 


r dr  = 


_ _ 1 n -2  1 (n+2)^ 


(5.2.14^ 


The  first  term  u T^  in  eq.  (5.2.14)  is  according  to  eq.  (B.63) 

2 


u T = 2y  u^  i f (R)  = - 2y  u^  1-  (n+2)  = - f 


2y  (n+2).  (5.2.15) 


Comparing  eqs.  (5.2.13)  and  (5.2.15) ,_w£  see  that  the  term  Wq  is 
indeed  small  relative  to  the  magnitude  of  u Tj.  In  Section  3 such  a 
ratio  of  magnitudes  was  anticipated  based  on  plausibility  arguments. 

The  total  correction  is  the  sum  of  Ca  and  Wq.  Combining  eqs. 
(5.2.13)  through  (5.2.15)  we  obtain  for  th 


e sum 


C.  + W = u f 
Co  z 


n-2  _ 2 

n+4 

f ^ y 3(n+lA 

2n  ~ \zj  3 

(n+1)  (n+2) 

^ y 4 (n+4) J 

(5.2.16) 


27 


The  right-hand  side  of  the  energy  equation  (5.2.12)  is  thus 


lT+<J=H-uf"  +C^  + W = H-  'uT  (1-a) 
z $ 0 z ^ 


(5.2.17) 


where  a is  a relative  correction  which  is  to  be  applied  to  u T^.  It  is 
given  by 

— fl  + (5.2.18) 

2)  p 4 n+4j  . 


The  second  term  in  this  formula  can  in  general  be  neglected,  because 
(R/z)^  is  of  the  ordor  lO”'^.  The  first  term  is  zero  only  for  n = 2, 
i.e.,  for  a Hagen-Poiseuille  flow  profile.  In  this  case  the  shear 
stress  is  a linear  function  of  r,  which  causes  certain  correction  terms 
to  vanish,  as  shown  in  Section  5.1.  For  a flat  flow  profile  with,  say, 
n = 10,  the  relative  correction  is  a = 0.4.  Clearly  such  a 40%  approx- 
imation error  will  be  seldom  tolerable.  Hence _for  f^a^  flow  profiles 
and  constant  viscosities  the  approximation  of  $by  -u  T2  is  not  realis- 
tic for  calculations  in  interior  ballistics. 

The  flow  profile  which  is  defined  by  eq.  (5.2.4)  does  not  have 
the  characteristic  form  of  a fully  developed  turbulent  flow  profile  for 
any  n.  We  may  therefore  ask  whether  the  correction  terms  are  possibly 
smaller  for  such  a profile.  In  order  to  investigate  this  question  we 
approximate  the  universal  turbulent  profile  (see,  e.g..  Reference  8, 
page  512)  by  defining 


n-2  /rV  2 n+4 

2n  ' 3 (n+l)(n+ 


f(r)  = 0.456 


The  corresponding  function  h(r)  is 


(5.2.19) 


h(r)  = - 0.456 


r 

R 


(5.2.20) 


Sahliahtingj  Boundary  Layer  TheovUt  MoGran-Hillt  ^ovk 
(4th  Edition) i 1960. 
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The  correction  term  of  the  momentum  equation  can  now  be 
computed  using  eq.  (5.2.5).  The  result  is 

= - 0.090  1^  (p  u2)  . (5.2.21) 


In  analogy  to  eq.  (5.2.8)  we  conclude  from  eq.  (5.2.21)  that  the 
momentum  flux  term  in  the  momentum  equation  should  be  increased  by  9% 
in  the  present  case. 

Assuming  as  before  constant  viscosities,  we  obtain  for  the  average 
heat  dissipation  function 


$ = 2ii  u ^ 
R 


^ 7.840  + f-2-]  (1.699  p + u’)_  (5.2.22) 


For  the  product  -u  Tj  we  obtain 


- u T = - 2p  u2  1=-  f (R)  = 2p  u2  • 14.364  . 


(5.2.23) 


The  right-hand  side  of  the  energy  equation  (5.2.12)  is  therefore 


H+(|)=H-uT 


1 - 0.454  + 


PI  0.059 


+ 0.588 


P' 


i 


(5.2.24) 


The_error  which  is  introduced  by  replacing  ^ by  -u  T-  is  about  45%  of 
|u  T^j.  As  in  the  previously  treated  case,  such  errors  will  be  seldom 
tolerable. 

We  may  conclude  from  these  examples  that  the  magnitudes  of  correction 
terms  are  essentially  the  same  for  flow  profiles  described  by  eq. (5.2.4) 
as  for  profiles  described  by  eqs.  (5.2.19)  and  (5.2.20).  Using  conven- 
tional tube  flow  equations,  e.g.  from  Reference  4,  for  interior  ballis- 
tics calculations,  one  introduces  errors  in  the  momentum  and  energy 
equations  which  are  of  the  order  of  9-50%  of  several  of  the  terms. 

The  examples  indicate  that  an  investigation  of  magnitudes  of  the 
correction  terms  is  necessary  whenever  average  flow  equations  are  used  . 
to  describe  non-steady  tube  flows. 
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6.  CONCLUSIONS 


Tube  flow  governing  equations  for  average  properties  differ  from 
one-dimensional  flow  equations.  The  differences  are  caused  by  the  fact 
that  averaging  of  functions  and  multiplication  of  functions  are  not 
commutative  operations.  The  magnitudes  of  the  differences  depend  on  the 
particular  problem.  If  the  unsteady  tube  flow  is  of  a type  which  is 
encountered  in  interior  ballistics,  then  several  terms  in  the  equations 
can  be  in  error  by  up  to  50%. 

One  consequence  of  the  various  correction  terms  in  the  equations  is 
that  the  continuity  and  momentum  equations  cannot  be  combined  to  yield 
a simple  equation  for  the  average  axial  velocity  component.  Instead, 
the  original  equation  for  the  average  axial  momentum  component  is  the 
simplest  form.  Correspondingly,  the  energy  equation  should  be  formulated 
for  the  internal  energy  per  unit  volume  instead  of  using  the  specific 
internal  energy. 

The  popular  approximation  of  the  heat  dissipation  function  by  the 
product  of  average  velocity  and  average  shear  stress  is  appropriate  only 
in  the  simplest  cases,  e.g.,  for  steady  flows  or  flows  with  a Hagen- 
Poiseuille  velocity  profile.  In  other  cases  the  approximation  can  be 
off  by  up  to  50%.  In  cases  of  more  complicated  flows  even  th's  sign  of 
the  approximation  can  be  wrong.  Hence  the  approximation  should  not  be 
used  unless  one  can  demonstrate  its  validity  in  the  particular  case  of 
application. 

Formulas  for  the  correction  terms  in  the  governing  equations  can 
be  derived  for  other  than  simple  tube  flows  following  the  outline  of 
this  paper.  The  derivations  which  are  presented  in  some  engineering  text- 
books neglect  important  first-order  terras.  The  apparent  success  of  the 
inaccurate  equations  for  the  treatment  of  tube  flows  is  probably  due  to 
the  fact  that  the  neglected  terms  are  small  or  vanish  for  steady  flows, 
for  which  most  comparisons  between  calculation  and  experiments  are  made. 
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APPENDIX  A 


FORMULAS  IN  CYLINDRICAL  COORDINATES 

In  Section  2 through  4 a convenient  cartesian  tensor  notation  was 
used  to  derive  all  formulas.  If  the  results  are  to  be  used  for  axially 
symmetric  tube  flows,  then  it  is  more  convenient  to  use  cylindrical 
coordinates.  In  this  appendix  we  express  the  important  formulas  in 
these  coordinates. 

Stokes  equation  for  the  stress  tensor  of  an  isotropic  fluid  can  be 
expressed  in  coordinate  independent  form  as  follows  (Reference  3,  page 
132;  Reference  2,  page  144) 


T = 2y  e + (m*  - -j  u)  div  u^ 


(A.l) 


where  x is  the  stress  tensor,  e is  the  strain  rate  tensor,  u*  is  the 
velocity  vector  of  the  fluid,  and  I is  the  unit  tensor.  The  viscosities 
u and  u'  in  eq.  (A.l)  need  not  be  constant,  i.e.,  the  fluid  under 
consideration  need  not  be  homogeneous.  However,  y as  well  as  y'  must 
be  positive  or  zero. 

Next  we  compute  the  work  rate  of  the  viscous  forces  acting  within 
an  arbitrary  volume.  To  this  end  we  compute  the  inner  product  of  the 
viscous  forces  V-x  with  the  velocity  vector  u*  and  integrate  over  the 
volume.  The  result  can  be  expressed  as  follows: 


/ (u*  • (Y'x))  dV  = ^ (u*  • (x'n))  dS  - /$dV 


(A.  2) 


In  eq.  (A. 2)  n is  a unit  vector,  orthogonal  to  the  surface  of  the  volume 
V and  pointing  inward,  and  ^ is  the  heat  dissipation  function  defined  by 


‘’'^race  = ‘"^‘trace  3 “>  “*  ''’trace  ' 


Because  div  u*  = expressed  as  follows: 


''""’trace  * ' 3 ''’mce  . 


(A.  4) 
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Eq.  (A. 4)  corresponds  to  eq.  (2.10)  in  cartesian  coordinates.  Eq. 

(A. 2)  corresponds  to  eq.  (2.13)  in  cases  where  the  body  forces  Xj 
are  absent. 

We  now  express  the  various  quantities  appearing  in  the  equations 
using  cylindrical  coordinates.  Let  the  coordinates  be  r,  (p  and  z. 
Components  of  vectors  and  tensors  we  denote  by  attaching  corresponding 
indexes  to  the  quantities.  Thus,  the  velocity  vector  u*  is 


u*  = (u^,  u^,  up 


(A.  5) 


The  strain  rate  tensor  e has  the  following  components  (Reference  2, 
Page  602) 


c 

rr 


3u  , 3u,  , 3u 

r _ <}>  1_  _ z 

3r’  ^<j)4>’"r34-  ^r  '^r’  ^zz  " 3z  ’ 


(A.  6) 


The  vector  V-t  has  the  components  Tr,  T^  and  T^,  representing  the  vis- 
cous forces  acting  in  the  three  coordinate  directions.  The  components 
are 


(V-t)^ 


r 3<fi 


> 


(V-T)^ 


IItJ.  + i—M  + + 2 - T 

3r  r 3 4)  3z  ^ r ' 


(Vt)^ 


rz 

3r 


r 8* 


zz 

3z 
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r 


rz 


(A.  7) 
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We  specialize  these  equations  for  the  case  of  an  axisymmetric  flow 
without  swirl  through  a circular  tube.  The  flow  is  then  independent 
of  the  coordinate  ({i,  and  the  (}> -component  of  the  velocity,  u^,  is 
zero.  In  order  to  simplify  the  notation  we  denote  the  non-zero 
velocity  components  as  fellows: 


= u(t,z,r) 


(A.  8) 


= v(t,z,r)  . 


Let  R be  the  radius  of  the  tube.  The  average  density  is  then  defined 
by 


2 ^ 

p(t,z)  = -T  / p(t,z,r)  r dr  . 
R o 


(A.  9) 


The  average  axial  velocity  is 


2 

u(t,z)  = —= / o(c,z,r)  I’(t,z,r)  r dr  . (A. 10) 

FA  ^ ^ 


R?y(t,z)  c 


The  local  continuity  equation  is 


|£..  (r  p V)  . |j  (p  u)  = 0 


CA.ll) 


The  corresponding  equation  for  the  averages  is 


I7  * {p  0)  = 0 . 


(A.12) 
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The  local  balance  of  momentum  is  expressed  by  the  following  two  differ- 
ential equations: 


+ i |-  (p  r v^)  + |-  (p  u v)  + = T 

at  r ar  •'  az  •’  ar  r 


a(p  u)  1 a . ^ a , 2.  ap  _ 

4^  * - (P  r u V)  . ^ (p  u ) * ^ 


(A. 13) 


The  third  momentum  equation  is  staisfied  identically  because  of  our 
symmetry  assumptions.  The  right-hand  sides  of  eqs.  (A. 13)  depend  on 
the  strain  rate  tensor  e by  eq.  (A. 7)  and  (A.l).  In  our  case  the  strain 
rate  tensor  has  the  following  components: 


3r  ’ ^ 


_ 1 ,, 

= — V , e = T—  . 

r zz  az  ' 


(A. 14) 


= ifll 

^rz  ~ 2 [_az  * ar 


The  divergence  of  the  velocity  vector  is 


div  u*  = + e, , + e,, 

trace  rr  99  zz 


= i i_  (r  V)  + ^ 
r ar  az  • 


(A. 15) 
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The  viscous  stress  tensor  t has  the  components 


T = 2ii 

rr 


3v 

3r 


(w' 


-f 


div  u* 
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rz 
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(p’ 


div  u*  , 


) (A. 16) 
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(fiZ 


= 0 
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zz 


= 2u 


3u 

3z 


(P' 


- y y)  div  u* 


The  right-hand  side  of  the  local  momentum  equations  (A. 13)  is 


, - 3t  , 

(V*t)  =T  = frx  ) + — — - T 

•'r  r r 3r  rr''  3z  r < 


CV-T)^  = T^=i|j(r 


} (A.  17) 
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Substituting  (A. 16)  into  (A. 17) , we  obtain 


Tr"Fl?(2ur|^^-  [( 


2 * If  (f  - I “J  ‘li'' 


n + i- 

3v  3u 

dZ 

^ 3z  3r 

1 3 r 3v  3u~l  3 r,  3U  , , 2 , 

= ^ ^ 37]  " ^ ^ - 3 


iv  u*J 


(A. 18) 


The  momentum  equation  for  the  averages  is 

3fp  u)  . 3 2,  3TS  ^ „ 

^ + ^(pu)+-rli-=T  + C , 
3t  3z  ^ 3z  z m 


(A. 19) 


where  the  average  pressure  p is  defined  by 

R 


■p  = — / p(t,z,r)  r dr 


0 


(A. 20) 


and  the  average  viscous  force  Tj  by 

R 


T = ” / T (t,z,r)  r dr 
^ R 0 ^ 


(A. 21) 


The  correction  term  C-,  in  eq.  (A.  19)  is 


■“m  3z 


The  equation  for  the  local  internal  energy  is 


■ 

R 

p u^ 

P 

u^  r dr\ 

R J 

. 

0 

- 

(A. 22) 


(r  p e V)  * (p  e u)  ♦ p [i|j.  (r  V)  . Ilj  . = Q . aiv  q 


(A. 23) 
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The  heat  dissipation  function '$ is  given  by  eq.  (A. 4).  Sub- 
stituting the  strain  rate  tensor  components  from  eq.  (A. 6)  into 
(A. 4),  we  obtain 


$ = 2li 

Cav\ 

2 1 

/av 

au\^  /vV 

|te) 

*2 

(az 

ar)  (r) 

+ (y*  - 

2 A 

3^^ 

f-- 

[r  8r 

(r  V) 

aJl^ 

* aj  * 

The  energy  equation  for  the  average  internal  energy  is 


a(p  e)  ^ 

at 


It  tP  5 0) . p 


+ C 


e2  ’ 


where 


2 ^ 

^ / p e r dr  , 
R 0 


H = -,  / (Q  - div  q)  r dr  , 
R 0 


and 


2 R 

$ = ^ $ r dr  , 

R 0 


dr 
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n 5) 


(p  u e)  r dr  . 


0 


(A. 24) 


(A. 25) 


(A.  26) 


(A.  2 7) 


(A. 28) 

(A.  29) 


(A.  30) 
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The  term -u  is  often  used  instead  of  $ in  eq.  (A. 25).  In  that  case 
the  equation  becomes 


3t 


3z 


u 


T + C , 
z el 


j.  r 


"e2 


+ C. 


(A.  31) 


The  additional  correction  terms  Wq  and  C$  are 

2 R 

W = -T  / div(T  u*)  r dr 
0 _2  ■' 

R 0 


(A.  32) 


and 


2 ^ 

C.  = u T =■  / u**(7*t)  r dr  . 

i z _2 

R o 


(A. 33) 


In  eq.  (A. 32)  we  have 


and 


div(T  u*) 


V + T 


Therefore 


r dr  = 


r dr  , (a. 34) 


where  div  u*  is  given  by  (A. 15). 
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Tlie  intogrand  in  (A. 33)  can  be  obtained  from  eq.  (A.17).  Carrying 
out  the  substitutions,  we  obtain  for  the  second  correction  term 


9 R 

^ ~ ~2  ^ ^ ^ = 
R 0 

Rr 


= UT  -^( 

^ rM 


- 3t 

13.  . , rz 

r 3r  W 3z 


- -t  1 - 

r ({)({._ 


+ u 


3t 


O 


13  zz 

7 3?  * — 


r dr  . 


(A.  35) 


The  separate  expressions  for  C,},  and  Wq  might  be  of  interest  for  the 
discussion  of  approximations.  Usually  C$  is  neglected  completely  and 
is  assumed  to  be  small  by  plausibilit;^  arguments.  The  total  correc- 
tion, which  is  caused  by  replacement  of  ^ by  -u  T^,  is  the  sum  of  Wq  and 
C^.  The  sum  is,  of  course. 


(A.  36) 


Eq.  (A. 36)  may  be  more  advantageous  for  actual  calculations  than  (A. 34) 
and  (A. 35)  because  it  does  not  contain  derivatives  of  the  viscosities. 
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APPENDIX  B 


LAGRANGE'S  APPROXIMATION  TO  INTERIOR  BALLISTICS  FLOW 

I^  is  plausible  to  assume  that  in  a gun  tube  the  average  axial  vel- 
ocity u(z,t)  of  each  cross-section  is  a linear  function  of  the  distance 
z from  the  breech  of  the  weapon.  Let  Zp(t)  be  the  location  of  the  pro- 
jectile and  Up(t)  = dzp/dt  be  its  velocity.  The  above-mentioned 
Lagrange's  approximation  is  then 


u(z,t)  = 


(B.l) 


In  the  classical  Lagrange's  approximation  (B.l)  is  supplemented  with 
the  assumption  that  the  gas  density  in  the  tube  is  a function  of  time 
only. 

In  this  appendix  we  shall  investigate  some  consequences  of  these 
assumptions.  Particularly  we  are  interested  in  finding  if  there  is  a 
three-dimensional  viscous  tube  flow  which  satisfies  Lagrange's  assumptions 

First  we  will  consider  flows  in  which  the  gas  density  is  a separable 
function  of  z,  t,  and  the  radial  coordinate  r: 


P(r,z,t)  = g(r)  • P(z)  • K(t) 


(B.2) 


Later  we  will  specialize  our  considerations  to  the  classical  Lagrange's 
approMir’^tion,  where  P(z)  and  g(r)  are  constants. 

We  assume  that  g(r)  is  non-dimensional  and  normalized  by 


~ / g(r)  r dr  = 1 . 


R^  0 


(B.3) 


The  product  of  the  other  two  functions  in  eq.  (B.2)  Ls  then  the  average 
density 


R 


p(z,t)  = P(z)  K(t)  - -j  I P(r,z,t)  r dr 

R 0 


(B.4) 
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In  eq.  (B.l)  the  variables  z and  t are  already  separated.  We  assume 
that  the  dependence  of  u on  r can  be  separated  also,  such  that 


u(r,z,t)  = fCr)  • u(z,t) 


(B.5) 


It  was  shown  in  Section  4 that  a reasonable  definition  of  the  aver- 
age axial  velocity  u in  terms  of  the  local  velocity  u is 


tr  = — • ^ / u pr  dr  . (B.6) 

0 

With  this  definition  of  u we  have  the  following  relation  between  the 
nondimensional  functions  f(r)  and  g(r): 

-y  / f(r)g(r)  r dr  = 1 . (B.7) 

- R 0 ' 


The  funtions  u(z,t)  and  p(z,t)  satisfy  the  continuity  equation  (4.6), 
i.e. , 


95-. 

at 


Q) 


0 . 


(B.8) 


Substituting  the  product  P*K  for  p into  eq.  (B.8)  and  the  expression 
(B.l)  for  u we  obtain 


nt)  • K'(t)  ^ K(t)  = 0 . 


(B.9) 


This  equation  has  solutions  of  the  form 


V = P 


0 


,m  z (o) 

z \ 


(m+1) 


z (t) 
P 


(B.IO) 


with  arbitrary  ra.  In  eq.  (B.IO)  p^  is  the  average  density  of  the  gas 
in  the  tube  at  time  t = 0.  For  m = 0 we  obtain  the  classical  Lagrange's 
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solution.  More  generally  we  may  assume  p to  be,  e.g.,  of  the  form 


P = 


A + -V  A , 
o m+1  ra  L 


m 


A + A I /'4.  \ I 
o m Iz^(t) 


z (o) 

■ R.: .. 
z (t) 
P 


(B.ll) 


with  arbitrary  m,  Aq  and  Ajjj.  For  physical  reasons  m ^ 0,  Aq  ^ 0 and 
Am  ^ "Aq* 

Next  we  investigate  the  radial  velocity  component  v(r,z,t).  The 
local  continuity  equation  is 

^ t (r  p V)  . - — (p  w)  * — (p  u)  = 0,  (8.12) 


where  w is  the  angular  velocity  component.  Let  w = 0 (no  swirll  and 


V = V(z,t)  • h(r)  . (B.13) 

Eq.  (B.12)  can  then  be  expressed  by 

• g(r)  + P V [r  g(r)  h(r)]  + g(r)  f(r)  [p  u]  = 0 . (B.14) 

Eliminating  3p/3t  from  eq.  (B.14)  with  the  aid  of  eq.  (B.8),  we  obtain 

[-l  + f(r)]  g(r)-g-|"^  + P V i ^[r  g(r)  h(r)]  = 0 . (^.15) 

This  equation  is  satisfied  by  the  functions 


and 


vU,t)  = RtH|^  (B.16) 

6 r . I 8 r ir  • (B-17) 

H r g ^ 
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Eqs.  (B.13),  (B.16),  and  (B.l?)  give  the  local  radial  velocity 
v(r,z,t,)  for  any  flow  profile  specified  by  p,  u,  g,  and  f.  Clearly, 
the  factor  h(r)  is  not  normalized  in  the  same  manner  a^  f(r).  There- 
fore, v(z,t)  is  not  an  "average"  radial  velocity.  If  p is  given  by 
eq.  (B.ll),  then 


7(z,t) 


(B.18) 


It  is  interesting  to  note  that  v is  not  zero  for  z = 0 and  z = Zp.  This 
is  an  indication  that  the  assumption  (B.5)  about  separation  of  variables 
for  the  axial  velocity  is  not  valid  in  the  vicinities  of  the  breech  and 
♦■he  projectile.  These  regions  we  will  therefore  exclude  from  our  consid- 
erations . 


In  summary  we  have  found  a flow  field  in  a cylindrical  tube  which 
satisfies  the  local  continuity  equation  and  Lagrange's  assumption  (B.l). 
The  flow  field  is  described  by  the  following  functions 


u 


z (t) 
P 


Up(t)  • f(r) 


p = p(z,t)  • g(r) 


> (B.19) 


v = V(z,t)  • h(r)  I 

If  one  specifies  u,  then  p is  given  by  eq.  (B.ll)  and  v is  given  by 
eq.  (B.18).  The  dependence  of  the  flow  field  on  r can  be  specified  by 
two  functions,  g(r)  and  f(r),  from  which  h(r)  is  then  computed  by  eq. 
(B.17).  The  function  g(r)  has  to  be  positive  for  0 ^r  ^ R and  normal- 
ized by  eq-  (B.3).  We  assume  also  that  g'(0)  = 0.  The  function  f(r) 
has  to  satisfy  the  conditions 


and 


f'(o)  = 0 

f(R)  = 0 . 


(B.20) 
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It  is  normalized  by  eq.  (8.7).  Hence  we  have  a total  of  three  condi- 
tions which  restrict  the  choice  of  f(r) . 

Instead  of  specifying  f(r)  we  may  also  specify  h(r).  The  function 
f(r)  is  then  given  by 

f(r)  = 1 - 7“^^  h(r)]‘  . (B.21) 


Tne  function  h(r)  has  to  satisfy  the  following  four  conditions 

h(0)  = 0 . 
h"(0)  = 0 , 
h(R)  = 0 , 
h'(R)  = 1/R  . 


The  flow  field  also  has  to  satisfy  the  momentum  equations.  The 
analysis  of  these  equations  is  more  complicated  because  it  involves,  in 
addition  to  the  velocity  and  density  functions,  the  pressure  function 
p(r,z,t)  and  the  viscosities  u and  y',  which  in  general  are  variable. 

We  have  tried  to  restrict  our  considerations  to  the  special  case  with 
constant  viscosities  and  classical  Lagrange's  approximation  (i.e., 
g(r)  =1).  We  have  found  that  the  flow  field,  as  defined  by  eq.  (B.19), 
does  not  satisfy  the  momentum  equations  in  this  special  case.  The  for- 
mulas for  correction  terms,  which  we  shall  derive  at  the  end  of  this 
appendix,  are  therefore  to  be  considered  as  approximations  only. 

If  g(r)  H 1,  then  the  flow  field  is  given  by 


u (t) 

u = Tr(z,t)f(r:  = z f(r)  , (B.23) 


Note  that  according  to  eq.  (B.24)  the  radial  velocity  component  v is 
independent  of  z for  the  classical  Lagrange's  approximation. 

The  divergence  of  the  velocity  u*  of  this  flow  is  (see  eq.  (A. 15)) 


div  u* 


i 3 Cry)  + 

r 3r  32 


3z  “ z (t)  * 
P 


(B.26) 


The  components  of  the  forces  caused  by  the  viscous  stress  tensor  are 
given  by  eq.  (A. 18).  In  our  case,  assuming  constant  viscosities,  we 
obtain  for  the  r-component 


T^  = 2y  V i [(r  h')'  . i h]  + p ||  f = 


3z 


[-2p  f . „ f]  11=  - w ^ f . 


(B.27) 


The  z-component  of  the  force  is 


T^  = iu  u [r  f']' 


_£i 


u z - [r  f']' 
P 


(B.28) 


The  local  momentum  equations  are  according  to  eqs.  (A. 13)  and  (B.12) 


3v  3v  ^ 3v  . 3p  ^ 

■rr  + py  -rr  + pu  TT  + -^r  " > 


^ Ji*  3r  • 3z  ■ 3r  *r 


(B.29) 


3u 


3u  ^ . 3u  ^ ^ T 
+ py  + pu  - T,  . 


3t  3r  3z  3z  z 


(B.30) 


For  a flow  field  described  by  eqs.  (B.23)  through  (B.24)  eq.  (B.29)  is 


pR 


t (--]  h+pR^  h h'  + 15-  = - u f . 


(B.31) 


Differentiating  eq.  (B.31)  with  respect  to  z,  we  obtain 


1^=0 

3r3z 


(B.32) 


The  function  p(r,z,t)  is  therefore  of  the  form 


p(r,z,t)  = Pj(r,t)  + . (B.33) 


Eq.  (B.31)  might  be  used  to  determine  the  function  p (r,t)  if  the  other 
terms  in  the  equation  are  given. 

Eq.  (B.30)  is  in  the  present  case,  i.e.,  for  the  flow  described  by 
eqs.  (B.23)  through  (B.25) 

or 

p VP/  p 

From  eq.  (B.35)  we  can  conclude  that  the  expressio:.  (3p2/'*z)/z  is  inde- 
pendent of  z.  The  various  terms  in  this  equation  are  products  of  func- 
tions of  r and  t and  the  equation  has  the  form 


g^(t)f^(r)  + g2(t)f2(r)  + gj(t)  + fj(r)  = 0 . (B.36) 


Such  an  equation  can  be  satisfied  identically  only  if  either  all  gi(t) 
are  constant  or  all  f^  (r)  are  constant.  The  case  with  all  fj^(r)  = const, 
corresponds  to  a slug  flow  in  which  we  are  not  interested.  Assuming  the 
time  functions  g^(t)  to  be  constant,  we  obtain  first 

goCt)  = p ^ = A 
^ P 
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or 


p 


Eq.  (B.37)  can  be  integrated  to  yield 


z (t)  = z (o)  ■ - ■ 

^ p"  1-At/p, 


The  corresponding  velocity  of  the  projectile  is 


u (t)  = z (o)  ^ r 

P P ^2 


0 (1-At/p^)' 


The  density  as  a function  of  time  is 


P(t)  =A^=p  fl-A^l  . 
P V 0. 


The  first  time  function  in  eq.  (B.35)  is  then 


Let  the  value  of  the  third  time  function  be  B.  We  obtain  then 


Eq.  (B.35)  takes  now  the  form 


A(f+f^+R  h f’)  + B - w - (r  f')»  = 


SO 


(B.37) 


(B.38) 


(B.39) 


(B.40) 


(B.41) 


(B.42) 


(B.43) 


The  functions  f(r)  and  h(r)  are  related  by  eq.  (B.21).  We  can  there- 
fore express  (B.43)  in  terms  of  h(r)  only.  We  also  multiply  the  equa- 
tion by  R^/y.  The  result  is 


with  the  constants 


and 


* 

A 


R u p 

2_  . 

y 


(B.45) 


(B.46) 


The  first  factor  in  eq.  (B.45)  is  a Reynolds  number  of  the  projec- 
tile. It  is  typically  of  the  order  10^.  The  first  factor  in  eq.  (B.46) 
for  B*  is  a Poiseuille  number  of  the  flow.  Its  magnitude  is  of  the 
order  10^.  The  function  h(r)  has  to  satisfy  the  differential  equation 
(B.44)  and  the  four  boundary  conditions  (B.22).  Since  eq.  (B.44)  is  of 
third  order  only,  the  function  h(r)  will  in  general  not  satisfy  all 
boundary  conditions.  We  conclude,  therefore,  that  Lagrange's  approxi- 
mation is  not  consistent  with  a flow  field  which  can  be  described  by 
separation  of  variables,  pqs.  (B.23)  through  (B.25). 


In  Section  5.2  we  have  nevertheless  used  this  flow  field  to  obtain 
estimates  of  correction  terms  because  we  were  not  able  to  find  an  exact 
three-dimensional  solution  of  Navier-Stokes  equations  which  is  also 
consistent  with  Lagrange's  approximation. 


Next  we  compute  the  various  correction  terms  for  the  average  flow 
equations  using  the  formulas  of  Appendix  A and  the  flow  described  by 
eq.  (B.19). 
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The  correction  terra  of  the  raoraentum  equation  is  given  by  eq. 
(A. 22): 


" 3z 


r 


L 


R 


-2  2 ( 
“ '7] 


pu  r dr 


J 


CB.47) 


Substituting  the  expressions  (B.23)  and  (B.25)  for  u and  p respectively 
into  eq.  (B.47),  we  obtain 


R 


'=m  “ It  [5  ^ '*’■] 


R 


-7i 


f^  r dri  1^  (p  u^) 


(B.48) 


Eq.  j[^B.48)  is  of  course  valid  for_  any  functions  u = u f(r)  and 

p = p(z,t).  In  Lagrange’s  case  p is  independent  of  z,  and  u is  linear 


C = 
m 


case 

r 

R 

1 “ h 

( f^  r dr 

R 

P 


(B.49) 


The  energy  equation  has  several  correction  terras.  First  we  consid- 
er the  term  given  by  eq.  (A. 29): 


We  substitute  the  flow  field  formulas  (B.19)  into  this  equation  and 
obtain 


R 


dr  = 


0 

R 
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This  expression  can  be  transformed  by  partial  integration  and  some 
algebra  into 

R R 

' f Hi  [?•-  h)'  P dr  - III]  (r  h)'  p dr  . CB.52) 


Eq.  (B.52)  shows  that  the  correction  term  Cei  is  zero  in  the  following 
cases : 

(a)  g = 1 and  |^  = 0,  i.e.,  the  classical  Lagrange's  case; 

o Z 
3p 

(b)  g = 1 and  — = 0,  i.e.,  p and  p independent  of  r; 

oT 

(c)  h = 0,  or  f = 1,  i.e.,  slug  flow. 

The  second  correction  term  of  the  energy  equation,  Cq2>  depends  on 
the  local  internal  energy.  According  to  eq.  (A. 30), 


■'e2  3z 


R 


2 

p u e : 

R^ 


j P u 


e r dr 


J 


(B.53) 


If  the  flow  field  is  given  by  eq.  (B.19),  then 

r R 


^e2  = i7^^ 


“"r  \ 

R P 3 


p e r dr  - 


R‘ 


3z 


L 

r 


L 


g (l-f)  e r dr 


J 


R 7 

^ puer  dr>  = 


J 


(B.54) 


This  correction  term  vanishes  if  the  internal  energy  e is  independent 
of  the  radial  coordinate  r. 


The  remaining  correction  terms.  Wo  and  C|j,  in  the  energy  elation 
are  caused  ^ the  replacement  of  the  heat  dissipation  function  by  the 
product  -u  Tj.  The  heat  dissipation  function  ^ is  according  to  eq. 

(A. 24) 


$ = 2ii 


2 2 
V^h' 


(B.55) 


53 


For  the  flow  field  described  by  eqs.  (B.23)  through  (B.25)  we  obtain 


$ 3 2y 


.£,)(\y.  CB.56) 

^ py 


The  average  (|>  is  by  definition 


$=  $r  dr  . 

R 0 


(B.57) 


The  correction  term  Wq  is  according  to  eq.  (A. 34) 

R,- 


° r2  J 


u 


c-  K i 1 > I-  (u’  - I P)  |S-E-  £ 

p ^ f V h r 2p  -jj-  u 3 ‘ 3zz_ 


r dr  = 


= 

Iv 


i 


h+2f^)  + (y'  - j p)  f r dr  . 


(B.58) 


The  correction  term  C,j,- is  given  by  eq.  (A, 35).  In  the  present  case 
with  constant  viscosities  we  obtain  by  substituting  eqs.  (B.27)  and 
(B.28)  into  eq.  (A. 35) 

R ^ 

-y^h-^f'+yufz-^  — (rf)  |rdr  = 


C,  " u T - \ ( 

‘ ^ R^J 


z r 

P P . 


0 

R,- 


= u T + 
z 


y R h f + y z^f’^ 


‘u^2 

^ r dr  . 

j v"p/ 


(B.59) 


By  partial  integration  and  using  eq.  (B.21)  we  can  show  that 


R 

/ R h f'r  dr 
0 


1 


r dr  . 


CB.60) 
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Making  use  of  eq.  (B.60)  we  can  express  the  correction  terms  as 
follows: 


0 


The  first  term  in  eq.  CB.62)  is  according  to  eq.  (B.28) 

u T = ^./  (r  f)'  dr  = 2ii  z \ f (R)  . (B.63) 

R 0 V.  py 

The  correction  terms_^  and  vanish  for  a slug  flow.  However, 
for  such  a flow  u T^  and  '$  are  also  zero. 
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Correction  term  in  energy  equation 
Definition  by  eq.  (4.23) 
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Definition  by  eq.  (4.24) 
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-1  -2 
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Radial  coordinate  (m) 
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Correction  terra  in  energy  equation  (Js"lm"2) 
Definition  by  eq.  (4.27) 


Cartesian  coordinate  (m) 


LIST  OF  SYMBOLS  (Cont'd) 

Strain  rate  tensor  (s“^).  Definition  by  eq.  (2.9) 
and  (A. 6) 

Ordinary  dynamic  viscosity  (.Pa*s) 

Dilatational  dynamic  viscosity  (Pa»s) 

Density  (kg/m^) 

Viscous  stress  tensor  (Pa) 

1 3 

Heat  dissipation  function  (J-s"'‘m“  ) 
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US  Array  Research  Office 
ATTN:  CRD-AA-EH 
P.  0.  Box  12211 
Research  Triangle  Park 
NC  27709 

1 Director 

US  Atihv  BMD  AGViiilc^Jci 

Technology  Center 
P.  0.  Box  1500,  West  Station 
Huntsville,  AL  35807 

1 Commander 

US  Army  Ballistic  Missile 
Defense  Systems  Command 
P.  0.  Box  1500 
Huntsville,  AL  35807 

3 Commander 

US  Naval  Air  Systems  Command 
ATTN:  AIR- 604 
Washington,  DC  20560 

3 Commander 

US  Naval  Ordnance  Systems 
ATTN:  ORD-9132 
Washington,  DC  20360 


1 Office  of  Naval  Research 
ATTN:  Code  460 
Washington,  DC  20360 

2 Commander  and  Director 
David  W.  Taylor  Naval  Ship 

Research  5 Development  Center 
ATTN:  Tech  Lib 

Aerodynamic  Lab 
Bethesda,  MD  20084 

4  Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Code  031,  Dr.  K.  Lobb 
Code  312,  Mr.  R.  Regan 
Mr.  S.  Hastings 
Code  730,  Tech  Lib 
Silver  Spring,  MD  20910 

4 Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Code  GX,  Dr.  W.  Kemper 
Mr.  F.  H.  Maine 

Ul  t VJ  » ^ 

Tech  Lib 

Dahlgren,  VA  22448 

3 Commander 

US  Naval  Weapons  Center 
ATTN:  Code  553,  Tech  Lib 

^ t . ^ 

UUUC 

Dr.  W.  Hazeltine 
Code  511,  Mr.  A.  Rice 
China  Lake,  CA  93555 

3 Commander 

US  Naval  Research  Laboratory 
ATTN:  Tech  Info  Div 

Code  7700,  D.  A.  Kolb 
Code  7720,  Dr.  E.  McClean 
Washington,  DC  20375 
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2 Commander  2 

US  Naval  Ordnance  Station 
ATTN:  Code  FS13A,  P.  Sewell 
Code  5235C,  A.  Horst 
Indian  Head,  MD  20640 

2 Commandant 

US  Marine  Corps 
ATTN:  A04F 

AX  1 

Washington,  DC  20380 

1 Commander 

Marine  Corps  Development  and 
Education  Center 

ATTN:  Chief,  Grounds  Oprs  Div  1 
Quantico,  VA  22134 

1 Marine  Corps  Liaison  Officer 
US  Army  Aviation  Test  Board 
Fort  Rucker,  AL  36362 

1 HQ  USAF/DPPTB  Maj  E.  G.  Nilson 
Washington,  DC  20330  1 

1 AFRPL/LKCB  (Dr,  Homing) 

Edwards  AFB,  CA  93523 

2 ADTC  (ADBPS-12) 

Eglin  AFB,  FL  32542 

3 AFATL  (DLDL)  1 

Eglin  AFB,  FL  32542 

2 AFATL  (DLDL,  Dr.  D.C.  Daniel; 

Tech  Lib) 

Eglin  AFB,  FL  32542 
1 AFWL  (DEV) 

Kirtland  AFB,  NM  87117  1 

1 ASD/XRA  (Stinfo) 

Wright-Patterson  AFB,  OH 
45433 


Director 

National  Aeronautics  and 
Space  Administration 
George  C.  Marshall  Space 
Flight  Center 
ATTN:  MS-1,  Lib 

R-AERO-AC,  Mr,  A.  Felix 
Huntsville,  AL  35812 

Director 

Jet  Propulsion  Laboratory 
ATTN:  Tech  Lib 
2800  Oak  Grove  Drive 
Pasadena,  CA  91103 

Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
ATTN:  MSI 85,  Tech  Lib 
Langley  Station 
Hampton,  VA  23365 

Director 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
ATTN:  MS  60-3,  Tech  Lib 
2100  Brookpark  Road 
Cleveland,  OH  44135 

Director 

NASA  Scientific  and  Technical 
Information  Facility 
ATTN:  SAK/DL 
P.  0.  Box  8757 
Baltimore/Washington 
International  Airport,  MD  21240 

Advanced  Technology  L aboratories 
ATTN:  Dr.  J.  Erdos 
Merrick  8 Stewart  Avenues 
Westbury,  NY  11590 
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2 ARO,  Inc. 

ATTN:  Tech  Lib 
Arnold  AFS,  TN  37389 

1 ARTEC  Associates.  Inc. 

ATTN:  Dr.  S.  Gill 
26046  Eden  Landing  Road 
Hayward,  CA  94545 

1 AVCO  Systems  Division 
ATTN:  Dr.  W.  Reinecke 
201  Lowell  Street 
Wilmington,  MA  01887 

1 Calspan  Corporation 

ATTN:  Mr.  G.  A.  Sterbutzel 
P.  0.  Box  235 
Buffalo,  NY  14221 

1 Technical  Director 

Colt  Firearms  Corporation 
150  Huyshore  Avenue 
Hartford,  CT  14061 

2 Dyna  East  Corporation 
ATTN:  Donald  L.  Tuckmantel 

Pei  Chi  Chou 
227  Hemlock  Road 
Wynnewood,  PA  19096 

1 General  Electric  Corporation 
Armaments  Division 
ATTN:  Mr.  R.  Whyte 
Lakeside  Avenue 
Burlington,  VT  05401 

1 Paul  Gough  Associates,  Inc. 
P.  0.  Box  1614 
Portsmouth,  NH  03801 

1 Martin  Marietta  Aerospace 
ATTN:  Mr.  A.  J.  Culotta 
P.  0.  Box  5387 
Orlando,  FL  32805 
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Copies  Organization 

1 Northrop  Corporation 
Aircraft  Division 
ATTN:  Dr.  A.  Wortman 
3901  W.  Broadway 
Hawthorne,  CA  90250 

1 Winchester-Western  Division 
Olin  Corporation 
New  Haven,  CT  06504 

1 Sandia  Laboratories 

ATTN:  Aerodynamics  Dept 

Org  5620,  R.  Maydew 
Albuquerque,  NM  87115 

1 S8D  Dynamics,  Inc. 

ATTN:  Dr.  M.  Soifer 
755  New  York  Avenue 
Huntington,  NY  11743 

1 Guggenheim  Aeronautical  Lab 
California  Inst  of  Technology 
ATTN:  Tech  Lib 
Pasadena,  CA  91104 

1 Drexel  University 

Dept  of  Mechanical  Engineering 
ATTN:  Prof.  Fa-Kang  Tsou 
Philadelphia,  PA  19104 

2 Franklin  Institute 
ATTN:  Dr.  Carfagno 

Dr.  Wachtell 
Race  8 20th  Streets 
Philadelphia,  PA  19103 

1 Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  MO  20810 
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1 Massachusetts  Institute  of 
Tecnnology 

Dept  of  Aeronautics  and 
Astronautics 
ATTN:  Tech  Lib 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 

1 Ohio  State  University 
Dept  of  Aeronautics  and 
Astronautical  Engineering 
ATTN : Tech  Lib 
Columbus,  OH  43210 

1 Penn  State  University 
A"TN:  Ken  Kuo 

York,  PA  17403 

2 Polytechnic  Institute  of 

Brooklyn 
Graduate  Center 
ATTN:  Tech  Lib 

Dr.  G.  Moretti 
Farmingdale,  NY  11735 

1 Princeton  University 
Forrestal  Campus  Library 
P.  0.  Box  710 
Princeton,  NJ  08540 


No.  of 

Copies  Organization 

1 Southwest  Research  Institute 
ATTN:  Mr.  Peter  S.  Westine 
P.  0.  Drawer  28510 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1 University  of  Illinois 

Aeronautical  and  Astronautical 
Engineering  Department 
ATTN:  H.  Krier 
Urbana,  IL  61801 

I University  of  Iowa 
College  of  Engineering 
ATiN:  Dr.  E.  Haug 
Iowa  City,  10  52242 

Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Dir,  USAMSM 
ATTN:  Dr,  J.  Sperrazza 
Mr.  Simmons 
Mr.  Clifford 
Dir,  USAMTD 
ATTN:  Small  Arms  Br 


1 Rensselaer  Polynic  Institute 
Dept  of  Mathematical  Sciences 
ATTN:  Dr.  Donald  Drew 
Troy,  NY  12181 
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